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Abstract 

We determine the 6-quark shape function parameters, A^^ and A^^ using the Behe B Xgj 
photon energy spectrum. We assume three models for the form of the shape function; Exponential, 
Gaussian and Roman. 
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INTRODUCTION 



The off-diagonal element Vub in the CKM matrix is extracted from measurements of the 
B — > Xuiu process in the limited region of lepton momentum or the hadronic recoil 
mass Mx 0] , or Mx and the leponic invariant mass squared where the contribution of 
background from the B — >■ Xciiy process is suppressed. In order to determine \Vub\ we need 
to extrapolate measured rates from such limited regions to the whole phase space. This 
extrapolation factor is evaluated using a theoretical prediction that takes into account the 
residual motion of the 6-quark inside the B meson, so called "Fermi motion" Q. Fermi 
motion is included in the heavy quark expansion by resumming an infinite set of leading- 
twist corrections into a shape function of the 6-quarkj^, 0]. Since the shape function is 
not calculable theoretically, it has to be determined experimentally. 

The best way is to make use of the photon energy spectrum for B Xg'j since both 
the inclusive decay spectra in _B ^ X„£z/ and B — > X^'y are expressed by the same shape 
function up to leading order of l/rrih in the heavy quark expansion The first results 

were obtained by CLEO but the errors of the shape function parameters are rather 
large. Therefore the uncertainty of the shape function dominates the theoretical error of 
\Vub\ at present. Belle has recently provided more precise data than CLEO of the B — > Xg'j 



photon spectrum We report on the results of determination of the shape function 



parameters using the Belle B — > Xgj data. 



PROCEDURE 

We used a method based on that devised by the CLEO collaboration 12]. We fit Monte 
Carlo (MC) simulated spectra to the raw data photon energy spectrum. "Raw" refers to the 
spectra that are obtained after the application of the B — »• Xgj analysis cuts. The use of 
"raw" spectra correctly accounts for the Lorentz boost from the B rest frame to the center 
of mass system, energy resolution effects and avoids unfolding. The method is as follows; 

1. Assume a shape function model. 

2. Simulate the photon energy spectrum for a certain set of parameters; (A^^, A^^). 

3. Perform a fit of the simulated spectrum to the data where only the normalization 
of the simulated spectrum is floated and keep the resultant value. 

4. Repeat steps 2-3 for different sets of parameters to construct a two dimensional grid 
with each point having a x^- 

5. Find the minimum on the grid and all points on the grid that are one unit of 
above the minimum. 

6. Repeat steps 1-5 for a different shape function model. 
Shape function models 

Three shape function forms suggested in the literature are employed; Exponential, Gaus- 
sian and Romanj^l^. These are described in Table B The shape function F is a function of 
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k+{= ko — fca), where fc^ is the residual momentum of the 6-quark in the B meson, defined 
through 

Pb,^l = rribV^ + k^, (1) 

where = (1,0,0,0) and k^ is the k component along the direction of the w-quark. The 
shape function is parameterized by A^^ and A^^. These parameters are related to the b 
quark mass, m^, and the average momentum squared of the b quark, /i^, via the relations, 

A^^ = MB-mb (2) 

and 

Af = -f^l (3) 

where Mb is the mass of the B meson. Up to leading order in the non-perturbative dynamics 
the shape function is universal in describing the 6-quark Fermi motion relevant to 6-to-light 
quark transitions. The lepton and photon energy spectra in i? — >■ and B —>■ Xs'j decays 
are given by the convolution of the respective parton-level spectra with the shape function. 
Example shape function curves are plotted in Figure Ufa). 



Shape Function 


Form 


Exponential 


F{k+;a) = iV(l - x)"e(i+")^ 


Gaussian 


F{k+; c) = N{1 - x)'=e-''(i-^)' 




where6=(F(^)/r(^))' 


Roman 


F{k+-p) = exp { - - k(1 - x)f} 




where k = -^eP/'^Ki{p/2) 




where x = k^/h?^ 




-mi, <k+< A^^ 




and a,c, p,N are chosen 




to satisfy 




Ao = 1, ^1 = 0, A2 = -Af /3, 




where An = J k\F{kj^)dk+ 



TABLE I: The three models used for the shape function forms 



Monte Carlo simulated photon energy spectrum 

We generate B —>■ Xg'j MC events according to the Kagan and Neubert prescription for 
each set of the shape function parameter values ^] . The generated events are then simulated 
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FIG. 1: Shape function model curves for Exponential {A^^,Xf) = (0.66,-0.40), Gaussian 
(A^^,Af^) = (0.63,-0.40), and Roman {A^^,\f) = (0.66,-0.39), where A^^" and are mea- 
sured in units of GeV/c^ and GeV^/c^ respectively. 

for the detector performance using the Belle detector simulation program. Afterwards B 
Xg'j analysis cuts are applied to the MC events to obtain the raw photon energy spectrum 
in the T(4S') rest frame jll|. 



Fitting the spectrum 

For a given set of shape function parameters, a fit of the MC simulated photon 
spectrum to the raw data spectrum is performed in the interval, 1.5 < E*/GeV < 2.8 [15j. 
The normalization parameter is floated in the fit. The raw spectrum is plotted in Figure |2l 
the errors include both statistical and systematic errors. The latter are dominated by the 
estimation of the BI3 background and are 100% correlated. Therefore the covariance matrix 
is constructed as 

■f/- ^stat^statx I ^sys^sys f a\ 

^ij = ^d, ^d^ ^ij + ^4 

where z, j = 1, 2, . . . , 13 denote the bin number, and cr^ is the error in the data. Then the 
used in the fitting is given by 

x^ = T.(d,-fm,r\d^-f,). (5) 

where {Vij)~^ denotes the ij^^ element of the inverted covariance matrix. The value after 
the fit is used to determine a map of as a function of the shape function parameters. 



The best fit and Ax^ contour 



The best fit parameters are associated to the minimum chi-squared case, Xmin- The la 
"ellipse" is defined as the contour which satisfies Ax^ = — Xmin = 1- The contours are 
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FIG. 2: Raw B Xg'j photon energy spectra in the T(45) frame as acquired from data. The 
errors include both statistical and systematic errors. Raw refers to spectra as measured after the 
application of Belle B Xs'j analysis cuts. 



found to be well approximated by the function [1 

The parameters a, b, c, d, and e are determined by fitting the function to the parameter 
points that lie on the contour. 



RESULTS 

The best fit parameters are given in Table |nj The parameter values are found to be 
consistent across all three shape function forms. The minimum fit for each shape function 
model is displayed in Figure |21 The fits to the contour with Ax^ = 1 points are shown in 
Figure^ The imposed shape function form acts to correlate A^^ and A^^. 



Shape 


A mm 


ASF 
(GeV/c2) 


Af 
(GeVVc^) 


Exponential 


4.883 


0.66 


-0.40 


Gaussian 


4.272 


0.63 


-0.33 


Roman 


5.020 


0.66 


-0.39 



TABLE IL The best fit shape function parameter values 
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''hope = Exp x7n* = 4.883/12 






(a) Exponential 



(b) Gaussian 



(c) Roman 



FIG. 3: The minimum ^ fits of MC simulated spectra to the raw data for each shape function 
model 
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FIG. 4: The fitted A^^ = 1 contours for each shape function model 
Strong Coupling 

The strong coupling constant, a^, is an input into the parton-level calculations for both 
B —>■ Xg'-f and B — > spectra. By default as{fi) is evaluated at the mass scale /i = m^. 
To investigate the systematic effect of this choice the analysis is redone for /i = mb/2 and 
II = 2mh in the case of the exponential shape function model. The A^^ and A^^ parameter 
values corresponding to Xmin given in Table IIIIl 



COMPARISON WITH CLEO 

The CLEO collaboration has provided points which lie on their equivalent Ax^ = 1 
contour for the case of an exponential shape function model^^l- The data points are slightly 
different from those given in the Gibbons' report since the present data now includes the 
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/ 1 

H' 


rv (ill 


^SF 


^1 






(GeV/c2) 


(GeVVc^) 


rrib 


0.210 


0.66 


-0.40 


mb/2 


0.257 


0.65 


-0.41 


2mi, 


0.177 


0.68 


-0.43 



TABLE III: The best fit parameters for various as using the exponential shape function model 



uncertainty in the BB background Monte Carlo normalization ^ . 

We fit the functional form given in equation EJ to their contour data and find excellent 
agreement (a = 2.378, h = -0.347, c = 0.178, d = -0.426, e = 0.256). The minimum 
point for the CLEO data corresponds to (A^^, A^^)exp = (0.545,-0.342). We compare the 
CLEO and Belle contours in Figure El The regions bounded by the contours marginally 
overlap. The uncertainty in the Belle result is much reduced with respect to that of CLEO. 

Unfortunately we can not produce a combined Ax^ = 1 contour of the two experiments 
since a precise map of Ax^ as a function of A^^ and A^^ is not currently available for CLEO. 
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FIG. 5: The fitted Ax^ = 1 contours for CLEO (blue) and Belle (red) assuming an exponential 
shape function form. 



SUMMARY 



We have determined the 6-quark shape function parameters, A^^ and Xf^, from fits of 
Monte Carlo simulated spectra to the raw Belle measured B — > Xg'j photon energy spectrum. 
Raw refers to the spectrum as measured after the application of analysis cuts. We used three 
models for the form of the shape function; Exponential, Gaussian and Roman. We found 
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the best fit parameters; (ASF,Af )exp = (0.66,-0.40), (A^^, A?^)Gauss = (0.63,-0.33), and 
(A^^, A^^)Roman = (0.66,-0.39), where A^^ and A^^ are measured in units of GeV/c^ and 
GeV^/c^ respectively. We also determined the Ax^ = 1 contours in the (A^^, A^^) parameter 
space for each of the assumed models. 
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